Methods and apparatus for calibrating a digital color imaging device that uses multi-hue colorants

ABSTRACT

Methods and apparatus are provided for calibrating a digital color imaging device to a printing press by determining a total colorant limit, per-channel colorant limits, and channel linearization tables using colorimetric and/or spectrophotometric techniques. In addition, for digital color imaging devices that use multi-hue colorants, methods and apparatus are provided for determining distribution functions for the multi-hue colorants as a function of input values.

CROSS-REFERENCES TO RELATED APPLICATIONS

This application is a continuation of U.S. patent application Ser. No. 10/836,069 filed Apr. 30, 2004, which is incorporated herein in its entirety by reference hereto.

BACKGROUND

This invention relates to color image processing. In particular, this invention relates to methods and apparatus for calibrating digital color imaging devices using calorimetric or spectrophotometric techniques. Digital color imaging devices, such as digital color printers and color copiers, have become increasingly popular in recent years. Indeed, while the cost of digital color imaging devices has dropped significantly, the number of hardware devices and software applications that are capable of producing color output that may be printed on such devices has substantially increased. Further, as the output quality and resolution of digital color imaging devices has improved, the number of uses for such devices has further increased.

For example, digital color laser printers and color inkjet printers now are increasingly used as relatively low cost proofing devices for commercial printing presses. Printing a print job on a printing press is a fairly expensive and time-consuming process. As a result, mistakes or errors in the print job are expensive to correct once a press run has commenced. To minimize such costly errors, high quality inkjet printers may be used to provide a proof of the print job before going to press. Ideally, the output of the proofing printer will visually match the output of the press. As a result, the proof output may be used for purposes of approving the print job or making any necessary modifications to the print job before printing the job on the press.

Referring now to FIG. 1, a previously known printing and proofing system is described. Printing system 20 includes commercial printing system 22 and proofing system 24. Commercial printing system 22 includes input device 26, input profile 28, color processing stage 30, press profile 32 and press 34. Input device 26 may be any device that may be used to create and/or store color image 38. For example, input device 26 may be a color scanner, digital camera, computer workstation, computer memory or other similar device.

Color image 38 includes a bitmap array of pixels, with each pixel including multiple colorant values. For example, if input device 26 is a scanner, color image 38 may include pixels expressed as a combination of red, green and blue (“RGB”) colorants. Colorant values typically are represented as multi-bit digital data values. Thus, if eight bits are used for each colorant, the colorant values may range from 0-255. In this regard, 0 corresponds to no colorant, and 255 corresponds to 100% colorant. The colorant values of color image 38 are defined in the device-dependent color space of input device 26.

Input profile 28 includes transformations between the color space of input device 26 and a profile connection space, such as Commission Internationale de I'Eclairage (“CIE”) XYZ, or other similar profile connection space. A profile connection space derived from the XYZ color space is commonly known as the CIE LAB color space, which expresses color values in a rectangular coordinate system, with the L, a, and b values each corresponding to one of the three dimensions in the system. The L-value characterizes the lightness aspect of the region along an axis ranging from black to white, with corresponding values ranging from 0 to 100. The a-value characterizes the color of the region along an axis ranging from green to red, with positive values corresponding to red and negative values corresponding to green. The b-value characterizes the color of the region along an axis ranging from blue to yellow, with positive values corresponding to yellow and negative values corresponding to blue. Together, the a-value and the b-value may be used to express the hue (“H”) and chroma (“C”) of the region:

$H = {\tan^{- 1}\left( \frac{b}{a} \right)}$ $C = \sqrt{a^{2} + b^{2}}$

The zero point in the plane defined by the a-values and the b-values corresponds to a neutral gray color having an L-value corresponding to the intersection of the plane with the L-axis.

Input profile 28 typically is produced in accordance with the profile specification of the International Color Consortium (“ICC”), and hence is referred to as an “ICC profile.” An ICC profile generally includes a transform from the profile connection space to the device space (the “forward transform”), and a transform from the device space to the profile connection space (the “backwards transform”). An input profile, however, typically includes only a backwards transform. For example, if input device is an RGB scanner, the backwards transform of input profile 28 may be used to convert device-dependent RGB colorant values to equivalent device-independent LAB colorant values.

Color processing stage 30 optionally may be used to perform various color processing operations in device-independent color space. For example, color processing stage may include software used to perform color editing or other color processing operations. Press profile 32 includes transformations between the color space of press 34 and a profile connection space, and also is typically an ICC profile. Thus, press profile 32 typically includes forward transform 32 a and backwards transform 32 b. For example, if press 34 is a conventional four-color offset press that uses cyan, yellow, magenta and black (“CMYK”) colorants, forward transform 32 a may be used to convert device-independent LAB colorant values to equivalent device-dependent colorant values CMYK₁. Press 34 receives CMYK₁ colorant values and provides press output 36 on media designed for use by a printing press.

Proofing system 24 includes press profile 32, printer profile 40, calibration stage 42 and proofing printer 46. In particular, backwards transform 32 b of press profile 32 may be used to convert device-dependent colorant values CMYK₁ to device-independent LAB colorant values. Printer profile 40 is typically an ICC profile, and includes a forward transform between the profile connection space and the color space of proofing printer 46. Accordingly, the forward transform of printer profile 40 is used to convert device-independent LAB colorant values to device-dependent colorant values CMYK₂.

Calibration stage 42 typically includes hardware and/or software that: (a) maps calibrated input values to equivalent uncalibrated input values (sometimes referred to as “linearization”); (b) limits the colorant of each channel; and (c) limits the total colorant of all channels. If proofing printer 46 uses multi-shade colorants, calibration stage 42 also may convert single colorant input values to equivalent multi-shade colorant values. The mapping and per-channel colorant limit functions typically are performed using tables that are designed to match the output response of proofing printer 46 to the output response of press 34, and also limit the colorant of each channel. The total colorant limit function is used to limit the total amount of colorant that may be output by proofing printer 46 to avoid negative image artifacts caused by using excessive colorant. Proofing printer 46 may be a digital inkjet printer, such as a CMYK inkjet printer or other similar printer. Proofing printer 46 receives calibrated CMYK colorant values and provides printed output 48 on media designed for use by an inkjet or laser printer.

The process of “calibrating” a printer typically includes determining linearization table values, per-channel colorant limits, a total colorant limit (“TCL”) and, optionally, distribution functions for multi-shade colorants. Referring now to FIG. 2, a previously known printer calibration process 50 is described. Beginning at step 52, a TCL is determined. In a multi-colorant printer, the amount of colorant for each channel typically is specified as a percentage between 0 and 100%. Thus, on a four-color printer, the maximum sum of all colorants that may be specified is 400%, corresponding to 100% on all four channels. If excessive colorant is used, however, undesirable image artifacts may result that produce an unacceptable print. For example, on inkjet printers, excessive colorant may cause bleeding (an undesirable mixing of colorants along a boundary between printed areas of different colorants), cockling (warping or deformation of the receiving material that may occur from using excessive colorant), flaking and smearing. In severe cases, excessive ink may cause the print media to warp so much that it interferes with the mechanical operation of the printer and may damage the printer. Thus, at step 52, a TCL is determined to minimize the effects of excessive colorant.

Previously known techniques for determining a TCL typically rely on trial and error methods that may be unsuitable for proofing purposes. In particular, previously known techniques typically involve printing several color patches that include various combinations of total amounts of colorant. A user then visually inspects the resulting printed output, and selects the patch (and thus the TCL) that produces the “best” results. A problem with such previously known techniques, however, is that the results may vary substantially from user to user, and even from time to time by the same user. The resulting lack of repeatability impairs the goal of obtaining a highly accurate proof.

Referring again to FIG. 2, after determining a TCL, at step 54 a colorant limit is determined for each channel. In a conventional printer, such as a CMYK inkjet or laser printer, the chroma response of the C, M and Y colorants as a function of the colorant amount is quasi-linear. However, beyond a certain specified colorant amount, the chroma actually begins to decrease, and the chroma response becomes highly non-linear. For the K channel, the luminance decreases with increasing colorant amount, until the luminance reaches a minimum level, but further increases in the colorant amount produce no further decrease in luminance. Indeed, for some combinations of colorants and media, oversaturation may occur, in which printed colors do not become any darker, and may actually become lighter, with increasing colorant amounts. Because it is difficult to accurately profile a printer in the non-linear region of operation, previously known techniques for calibrating a printer typically limit the colorant of each channel so that the printer operates only in the quasi-linear region and not in the oversaturation region.

Previously known techniques for determining per-channel colorant limits, however, have typically relied on density-based measurements that may be incomplete and inaccurate for proofing purposes. In particular, previously known techniques for determining per-channel colorant limits typically involve printing a target for each colorant, where the target includes several color patches that range from 0 to 100% colorant. After printing the target, a user typically measures the optical density of each patch using a densitometer or other instrument that provides optical density values. The per-channel colorant limits are then specified as the colorant values that produce a predetermined density (e.g., the lowest maximum density) on all channels.

One problem with such isometric density techniques is that they fail to consider the impact of the colorant limitation on the gray balance of the printer. When a printer outputs approximately equal percentages of C, M and Y colorants, a neutral gray should result. The human eye is very sensitive to detecting shifts in neutrality when neutral areas are compared side-by-side. Thus, gray balance may be used to determine if the gamut of one printing device (e.g., a proofing printer) matches the gamut of another printing device (e.g., a press). Previously known density-based techniques for determining per-channel colorant limits, however, typically do not ensure proper gray balance. To solve this problem, experienced users have developed their own techniques for achieving a desired density value for each colorant and also a good gray balance. Such empirical techniques vary from user to user, however, and require specialized knowledge that all users may not possess.

In addition, previously known density-based techniques for determining per-channel colorant limits may be inaccurate for proofing printers. Conventional densitometers typically operate by illuminating a printed patch using light having a known spectral distribution, and then measuring the amount of light absorbed in a narrow frequency band of the visual spectrum. The absorption measurement may then be translated to a density measurement, with higher absorption corresponding to higher density. Densitometers typically use narrow-band optical filters that are tailored to describe the behavior of colorants used on a conventional printing press. Unfortunately, however, the filters are not optimized for describing the behavior of colorants used by conventional inkjet and laser printers used for proofing. Indeed, if a colorant used by a proofing printer has a maximum absorption at a frequency outside the band of the instrument's filters, the resulting density measurements may be incorrect. As a result, density-based techniques for determining per-channel colorant limits may produce inaccurate results.

Referring again to FIG. 2, after per-channel colorant limits have been determined, in step 56, linearization tables are calculated for each channel so that the output response of proofing printer 46 matches the output response of press 34. Previously known techniques for calculating linearization tables typically involve printing a target for each colorant that includes several color patches that range from 0 to 100% colorant coverage. After printing the target, a user typically measures the optical density of each patch using a densitometer or other instrument that provides optical density values, and then calculates table values that map the input/output density response of the printer to an input/output density response of the press. As described above, however, conventional densitometers and similar measuring instruments may not accurately measure density of colorants used by conventional inkjet and laser printers used for proofing. As a result, previously-known density-based techniques for calculating linearization tables may produce similarly inaccurate results.

Referring again to FIG. 2, after linearization tables have been calculated, at optional step 58, distribution functions may be determined for multi-hue colorants. In particular, high-quality digital inkjet printers used for proofing purposes often include four primary CMYK colorants (also referred to herein as “normal cyan,” “normal magenta,” “normal yellow” and “normal black”), plus light cyan and light magenta colorants (indicated by lowercase “c” and “m”) to provide improved image quality in the highlight regions of an image. Referring again to FIG. 1, if proofing printer 46 is a CcMmYK printer, printer profile 40 typically converts LAB values to CMYK values, and calibration stage 42 converts cyan values into mixtures of normal cyan (C) and light cyan (c) values, and converts the magenta values into mixtures of normal magenta (M) and light magenta (m) values.

Previously known techniques for converting a specified colorant value to equivalent multi-shade colorants often rely on trial and error techniques to determine the distribution function between the colorants. The resulting distribution function may be acceptable for a first set of colorants (e.g., used in a first printer in location A), but may be unacceptable for a second set of colorants (e.g., use in a second printer in location B). As a result, unless a new distribution function is determined for the second set of colorants, the printed output of the two printers may not match. Previously known trial and error techniques, however, typically do not permit easy modification of distribution functions. Instead, the entire process must be repeated, which may be extremely time consuming and inefficient.

In view of the foregoing, it would be desirable to provide apparatus and methods for calibrating a digital imaging device in a repeatable manner.

It also would be desirable to provide apparatus and methods for calibrating a digital imaging device in an accurate manner.

It additionally would be desirable to provide apparatus and methods for calibrating a printer without requiring specialized knowledge by a user.

SUMMARY

In view of the foregoing, it is an object of this invention to provide apparatus and methods for calibrating a digital imaging device in a repeatable manner.

It also is an object of this invention to provide apparatus and methods for calibrating a digital imaging device in an accurate manner.

It additionally is an object of this invention to provide apparatus and methods for calibrating a printer without requiring specialized knowledge by a user.

These and other objects of this invention are accomplished by providing methods and apparatus for calibrating a digital color imaging device to a printing press by determining a total colorant limit, per-channel colorant limits, and channel linearization tables using calorimetric and/or spectrophotometric techniques. In addition, for digital color imaging devices that use multi-hue colorants, methods and apparatus of this invention optionally may determine distribution functions for the multi-hue colorants as a function of input values.

BRIEF DESCRIPTION OF THE DRAWINGS

The above-mentioned objects and features of the present invention can be more clearly understood from the following detailed description considered in conjunction with the following drawings, in which the same reference numerals denote the same elements throughout, and in which:

FIG. 1 is a block diagram of a previously known printing system;

FIG. 2 is a flowchart of a previously known printer calibration process;

FIG. 3 is a block diagram of an exemplary calibration system in accordance with this invention;

FIG. 4 is a flowchart of an exemplary process for determining a total colorant limit in accordance with this invention;

FIG. 5 is an exemplary test chart for use with the method of FIG. 4;

FIGS. 6A-6C are tables of exemplary colorant values and associated colorimetric measurements for the test chart of FIG. 6;

FIG. 7 is a flowchart of an exemplary process for determining per-channel colorant limits in accordance with this invention;

FIGS. 8A-8C are flow diagrams of exemplary methods for determining minimum, maximum and optimal per-channel colorant limits in accordance with this invention;

FIG. 9 is an exemplary test chart for use with the methods of FIG. 8;

FIGS. 10A-10G are tables of exemplary colorant values and associated calorimetric measurements for the test chart of FIG. 9;

FIGS. 10H-10I are tables of exemplary interpolated colorant values and calculated chroma values in accordance with this invention;

FIG. 11 is a flowchart of an exemplary process for determining a linearization table in accordance with this invention;

FIG. 12 is an exemplary test chart for use with the method of FIG. 11;

FIG. 13 is a tables of exemplary colorant values and associated calorimetric measurements for the test chart of FIG. 12;

FIG. 14 is a diagram of exemplary tonal responses calculated in accordance with the method of FIG. 11; and

FIG. 15 is a flow diagram of an exemplary process for determining dual-tone distributions in accordance with this invention.

DETAILED DESCRIPTION

Referring now to FIG. 3, an exemplary system in accordance with this invention is described for calibrating a digital imaging device using calorimetric and/or spectrophotometric techniques. Calibration system 70 includes image source 72, proofing printer 46, output pages 74, measurement device 76, press profile 32 and processor 78. Image source 72 includes image file 80 that comprises digital data representing test patterns 82 to be printed by proofing printer 46. Image source 72 may be a personal computer, laptop computer, handheld computer, computer workstation, print server, personal digital assistant, or any other similar device that may be used to provide image files for printing by color imaging devices.

Image source 72 may include a software application (not shown) used to generate image file 80. For example, image source 72 may be a personal computer that includes software that may be used to generate image file 80. Image file 80 may be a digital data file that describes test patterns 82 in a page description language, such as PostScript, PCL, or other similar page description language, or may simply be a raster image, such as a TIFF image, RAW image, or other similar raster image. Proofing printer 46 may be a laser printer, inkjet printer or other similar color imaging device that uses one or more colorants to provide output pages 74 including test patterns 82. For example, proofing printer 46 may be a four-color inkjet printer that uses CMYK colorants, a six-color inkjet printer that uses CcMmYK, or other similar multi-colorant imaging device. Test patterns 82 include one or more color patches P.

Measurement device 76 may be any conventional measurement device that may be used to provide spectral and/or colorimetric data that describes a printed sample, such as a colorimeter, spectrophotometer, spectrocolorimeter, or other similar device. For example, measurement device 76 may be a Spectrolino spectrophotometer manufactured by GretagMacbeth LLC, New Windsor, N.Y., U.S.A. Measurement device 76 provides colorimetric data, such as CIE LAB data (referred to herein as “LAB data”), CIE XYZ data (referred to herein as “XYZ data”), CIE LCH data (referred to herein as “LCH data,” where L-values correspond to lightness, C-values correspond to chroma, and H-values correspond to hue), or other similar colorimetric and/or spectral data that describes printed samples, such as color patches P.

Processor 78 may be a personal computer, laptop computer, handheld computer, computer workstation, print server, personal digital assistant, or any other similar device that may be used to receive colorimetric data, such as LAB data (i.e., L-, a- and b-values), LCH data (i.e., L-, C- and H-values), or other similar calorimetric and/or spectral data from measurement device 76 and generate therefrom calibration data. Persons of ordinary skill in the art will understand that the functions of processor 78 may be implemented by image source 72. In accordance with this invention, processor 78 determines total colorant limit, per-channel colorant limits, and channel linearization tables (one table per colorant) using calorimetric and/or spectrophotometric techniques. For proofing printers that use multi-hue colorants (e.g., light and normal cyan), processor 78 may also determine distribution functions for the multi-hue colorants as a function of input values. Each of these techniques will be discussed in turn.

Total Colorant Limit

Referring now to FIGS. 3 and 4, exemplary methods and apparatus in accordance with this invention are described for determining a TCL of proofing printer 46. Beginning at step 90, proofing printer 46 is used to print test pattern 82 including test patches P on output page 74. For example, a user of image source 72 may issue a print command to print image file 80 including test pattern 82 on proofing printer 46. An exemplary output page 74 a including exemplary test pattern 82 a is illustrated in FIG. 5. Test pattern 82 a includes an array of three strips A, B and C of test patches P, with each strip including sixteen test patches. Persons of ordinary skill in the art will understand that test pattern 82 a may include more or less than three strips, and each strip may include more or less than sixteen test patches P. Each test patch P is comprised of a corresponding specified percentage of colorants used by proofing printer 46 (e.g., C, M, Y and K).

FIGS. 6A-6C illustrate tables of exemplary colorant values (in percent) for test patches P in strips A, B and C, respectively (patches are identified in each table by row (A, B or C) and column number (1-16)). The exemplary values provide 48 test patches P having total colorant values (i.e., the sum of the colorant percentages for each patch) that range from 99% to 300% for strip A, 100% to 200% for strip B, and 100% to 400% for strip C. In this regard, the 48 exemplary test patches P cover a broad range of total colorant values between 99% and 400%. Persons of ordinary skill in the art will understand that other colorant values also may be used for test patches P, and that the range of colorant values may include values less than 100%. In addition, persons of ordinary skill in the art will understand that the arrangement of patches within strips A, B and C, and the arrangement of the strips within test pattern 82 a may be changed.

Referring again to FIGS. 3 and 4, at step 92, calorimetric values are determined for each test patch P printed in step 90. For example, measurement device 76 may be used to determine LAB data for each test patch P on output page 74 a. FIGS. 6A-6C illustrate exemplary measured LAB data for test patches P in test pattern 82 a. Referring again to FIG. 4, at step 94, a test patch P is identified in each strip A, B and C that has the minimum L-value of all of the patches in the strip. Thus, from the exemplary calorimetric values shown in FIGS. 6A-6C, the minimum L-value for strip A (10.53) corresponds to test patch A11, the minimum L-value for strip B (26.69) corresponds to test patch B9, and the minimum L-value for strip C (9.18) corresponds to test patch C9.

Referring again to FIG. 4, at step 96, the total area coverage (“TAC”) is determined for each of the patches identified at step 94. The TAC of a patch equals the sum of the colorant values for the patch. Thus, referring to FIGS. 6A-6C, the TAC for exemplary test patches A11, B9, and C9 is 300%, 200% and 260%, respectively. Referring again to FIG. 4, at step 98, the TCL is set to the maximum TAC determined in step 96. Thus, for the exemplary colorant values shown in FIGS. 6A-6C, the TCL is 300%.

Per-Channel Colorant Limit

In addition to determining TCL, methods and apparatus in accordance with this invention also determine a limit for each colorant of proofing printer 46, while seeking to maintain the gamut of the proofing printer as large as the gamut of press 34. Referring now to FIGS. 3 and 7, exemplary methods and apparatus in accordance with this invention are described for determining per-channel colorant limits. Beginning at step 100, a minimum limit is determined for each colorant. Next, at step 102, a maximum limit is determined for each colorant. Finally, at step 104, an optimal limit between the minimum and maximum limit is determined for each colorant. Each of these steps will be described in turn.

Referring now to FIGS. 3 and 8A, an exemplary method 100 for determining a minimum per-channel colorant limit is described. In particular, at step 110, proofing printer 46 is used to print test pattern 82 including test patches P on output page 74. An exemplary output page 74 b including exemplary test pattern 82 b is illustrated in FIG. 9. Test pattern 82 b includes an array of eight strips A-H of test patches P, with each strip including twenty-two test patches. Persons of ordinary skill in the art will understand that test pattern 82 b may include more or less than eight strips, and each strip may include more or less than twenty-two test patches P. Each test patch P is comprised of a corresponding specified percentage of colorants used by proofing printer 46 (e.g., C, M, Y and K).

FIGS. 10A-10G illustrate exemplary colorant values (in percent) for test patches P (patches are identified in each table by row (A-H) and column number (1-22)). The exemplary values provide test patches P having scales of single-colorant values for each colorant, and patches having various combinations of multiple-colorant values. In particular, FIG. 10A illustrates exemplary colorant values for test patches A1-A17 including a scale of cyan colorant from 100% to 28%; FIG. 10B illustrates exemplary colorant values for test patches A18-B10 including a scale of magenta colorant from 100% to 36%; FIG. 10C illustrates exemplary colorant values for test patches B11-C2 including a scale of yellow colorant from 100% to 40%; FIG. 10D illustrates exemplary colorant values for test patches H1-H22 including a scale of black colorant from 36% to 100%; and FIGS. 10E-10G illustrate exemplary colorant values for test patches C3-G22 including various combinations of C, M and Y colorants. Persons of ordinary skill in the art will understand that other specific colorant values also may be used for test patches P.

Referring again to FIG. 8A, at step 112, calorimetric values are determined for each test patch P printed in step 110. For example, measurement device 76 may be used to determine XYZ and LAB data for each test patch P on output page 74 b. FIGS. 10A-10G illustrate exemplary measured XYZ and LAB data (and corresponding chroma data) for test patches P. Referring again to FIG. 8A, at step 114, colorimetric values associated with maximum amounts of each colorant are extracted from press profile 32. In particular, for cyan, magenta, yellow and black, the colorimetric values that correspond to maximum amounts of cyan (i.e., CMYK value of (100, 0, 0, 0)), magenta (i.e., CMYK value of (0, 100, 0, 0)), yellow (i.e., CMYK value of (0, 0, 100, 0)) and black (i.e., CMYK value of (0, 0, 0, 100)), respectively, are extracted from backwards transform 32 a of press profile 32. Table 1 illustrates exemplary colorimetric values extracted from exemplary press profile 32:

TABLE 1 CMYK L a b (100, 0, 0, 0) 54 −36 −50 (0, 100, 0, 0) 47 75 −7 (0, 0, 100, 0) 88 −6 95 (0, 0, 0, 100) 18 1 −1

Referring again to FIG. 8A, at steps 116 a-d, the minimum limits C_(MIN), M_(MIN), Y_(MIN) and K_(MIN) are determined for the C, M, Y and K colorants, respectively. In particular, at step 116 a, the minimum cyan limit C_(MIN) is identified as the colorant value of the minimum cyan patch that has a- and b-values whose magnitudes are greater than or equal to the magnitudes of the a- and b-values, respectively, extracted for cyan from press profile 32. At step 116 b, the minimum magenta limit M_(MIN) is identified as the colorant value of the minimum magenta patch that has: (a) a- and b-values whose magnitudes are greater than or equal to the magnitudes of the a- and b-values, respectively, extracted for magenta from press profile 32, or (b) an L-value that is less than the L-value extracted for magenta from press profile 32, whichever is lower. At step 116 c, the minimum yellow limit Y_(MIN) is identified as the colorant value of the minimum yellow patch that has a b-value whose magnitude is greater than or equal to the magnitude of the b-value extracted for yellow from press profile 32. At step 116 d, the minimum black limit K_(MIN) is identified as the colorant value of the minimum black patch that has an L-value that is less than the L-value extracted for black from press profile 32.

For example, referring to FIG. 10A, cyan test patch A12 is identified as the minimum cyan patch that has calorimetric values a=−39.72 and b=−50.34 whose magnitudes are greater than the corresponding colorimetric values extracted from press profile 32 for maximum cyan (i.e., a=−36 and b=−50 (Table 1)). The cyan value of test patch A12 is 47%. Therefore, C_(MIN)=47%. Referring to FIG. 10B, magenta test patch B6 is the minimum magenta patch that has calorimetric values a=75.23 and b=−21.83 whose magnitudes are greater than the magnitudes of the corresponding calorimetric values extracted from press profile 32 for maximum magenta (i.e., a=75 and b=−7 (Table 1)). In addition, magenta test patch B1 is the minimum magenta patch that has calorimetric value L=46.58 that is less than the corresponding calorimetric value extracted from press profile 32 for maximum magenta (i.e., L=47). The magenta value of test patch B1 is 75% and B6 is 53%. Therefore, the lower of these two values is M_(MIN)=53%. Referring to FIG. 10C, yellow test patch B15 is the minimum yellow patch that has a calorimetric value b=95.89 whose magnitude is greater than the magnitude of the corresponding calorimetric value extracted from press profile 32 for maximum yellow (i.e., b=95 (Table 1)). The yellow value of test patch B15 is 80%. Therefore, Y_(MIN)=80%. Referring to FIG. 10D, black test patch H3 is the minimum black patch that has a calorimetric value L=16.80 that is less than the corresponding calorimetric value extracted from press profile 32 for maximum black (i.e., L=18 (Table 1)). The black value of test patch H3 is 42%. Therefore, K_(MIN)=42%.

Referring now to FIG. 8B, an exemplary method 102 for determining a maximum per-channel colorant limit is described. Steps 110 and 112 have already been described above in connections with FIG. 8A. At steps 118 a-d, the maximum limits C_(MAX), M_(MAX), Y_(MAX) and K_(MAX) are determined for the C, M, Y and K colorants, respectively. For the C, M and Y channels, the maximum limit C_(MAX), M_(MAX) and Y_(MAX), is identified as the colorant value of the minimum cyan, magenta and yellow patch, respectively, where: (a) the chroma is maximum, or (b) oversaturation begins. The maximum limit K_(MAX) is identified as the colorant value of the minimum black patch where: (a) the L-value is minimum, and (b) the difference between the L-value of the patch and the L-value of the 100% patch (L₁₀₀) is sufficiently small. For example, the difference is sufficiently small if the following inequality is satisfied:

(Int[L(i)]−Int[L ₁₀₀])<1  (1)

where Int[L(i)] is the integer portion of L(i) and Int[L₁₀₀] is the integer portion of L₁₀₀.

For C, M and Y colorants, chroma increases with increasing colorant value, but may begin to decrease beyond a certain colorant amount. As a result of rounding errors and device tolerances in measurement equipment, measured chroma values may fluctuate slightly about a maximum value. To avoid selecting a false maximum, therefore, the chroma value of each test patch is compared to the chroma values of several successive patches. For example, for a series of chroma values, the maximum chroma value may be identified as the chroma value that satisfies at least one of the following inequalities:

chroma(i)>chroma(i+1)>chroma(i+2)  (2a)

or

chroma(i)>chroma(i+1)>chroma(i+3)  (2b)

In this case, the maximum chroma value must be greater than the first and second successive chroma values, or the first and third successive chroma values.

To determine if C, M or Y colorants are oversaturated, X, Y and Z colorimetric values may be used. In particular, because X-values are proportional to red reflection, and because cyan absorbs red, X-values ideally decrease with increasing amounts of cyan colorant. If cyan becomes oversaturated, however, X-values may begin to increase with increasing cyan colorant. As a result, the start of oversaturation of cyan may be identified by determining the colorant value of the patch where X-values begin to increase with increasing cyan colorant. Similarly, Y-values ideally decrease with increasing amounts of magenta colorant, and Z-values ideally decrease with increasing amounts of yellow colorant. As a result, the start of oversaturation of magenta may be identified by determining the colorant value of the patch where Y-values begin to increase with increasing magenta colorant, and the start of oversaturation of yellow may be identified by determining the colorant value of the patch where Z values begin to increase with increasing yellow colorant.

As a result of rounding errors and device tolerances in measurement equipment, X, Y and Z measurements may fluctuate slightly about a minimum value. To avoid selecting a false minimum, therefore, the X-value of each test patch is compared to the X-values of several successive patches. For example, the start of oversaturation of cyan may be identified as the colorant patch having an X-value X(i) that satisfies at least one of the following inequalities:

X(i)<X(i+1)<X(i+2)  (3a)

or

X(i)<X(i+1)<X(i+3)  (3b)

In this case, the X-value must be less than the first and second successive X-values, or the first and third successive X-values. Similarly, the start of oversaturation of magenta may be identified as the colorant patch having an Y-value Y(i) that satisfies at least one of the following inequalities:

Y(i)<Y(i+1)<Y(i+2)  (4a)

or

Y(i)<Y(i+1)<Y(i+3)  (4b)

and the start of oversaturation of yellow may be identified as the colorant patch having a Z-value Z(i) that satisfies at least one of the following inequalities:

Z(i)<Z(i+1)<Z(i+2)  (5a)

or

Z(i)<Z(i+1)<Z(i+3)  (5b)

For K, the minimum L-value may be identified as the L value L(i) that satisfies the equation:

L(i)≦L(i+1)  (6)

Referring again to FIGS. 10A-10D, the above principles may be used to determine C_(MAX), M_(MAX), Y_(MAX) and K_(MAX). In particular, referring to FIG. 10A, the chroma value of cyan test patch A8 satisfies equation (2a), and therefore the maximum chroma (68.49) corresponds to a cyan colorant value of 65%. In addition, none of the X-values satisfy equation (3), and therefore cyan is not oversaturated. As a result, C_(MAX)=65%. Referring to FIG. 10B, the chroma value of magenta test patch A22 satisfies equation (2a), and therefore the maximum chroma (82.55) corresponds to a magenta colorant value of 80%. In addition, none of the Y-values satisfy equation (4), and therefore magenta is not oversaturated. As a result, M_(MAX)=80%. Referring to FIG. 10C, the chroma value of all yellow test patches satisfy equation (2a), and therefore the maximum chroma (102.15) corresponds to a yellow colorant value of 100%. In addition, none of the Z-values satisfy equation (5), and therefore yellow is not oversaturated. As a result, Y_(MAX)=100%. Referring to FIG. 10D, the L-value of black test patch H12 satisfies equations (1) and (6), and therefore the maximum limit for black is K_(MAX)=69%. Thus, following are minimum and maximum limits from the exemplary colorant data of FIGS. 10A-10D:

Colorant Min Limit Max Limit Cyan 47% 65% Magenta 53% 80% Yellow 80% 100%  Black 42% 69%

Referring now to FIG. 8C, an exemplary method 104 for determining an optimal per-channel colorant limit is described. Steps 110 and 112 have already been described above in connections with FIG. 8A. Beginning at step 120, the optimal limit C_(OPT), M_(OPT), and Y_(OPT) is determined for the C, M and Y colorants, respectively. In particular, colorant values are determined for a plurality of “interpolated patches,” consisting of combinations of the cyan, magenta and yellow colorants used to determine the minimum and maximum colorant values determined in steps 100 and 102. Thus, in the examples illustrated in FIGS. 10A-10C, interpolated patches are created using combinations of the following cyan, magenta and yellow colorants:

Cyan: 47%, 53%, 56%, 60% and 65% Magenta: 53%, 56%, 60%, 65%, 70%, 75% and 80% Yellow: 80%, 85%, 90%, 95% and 100%

Colorant values for such interpolated patches are illustrated in FIGS. 10H-10I.

Referring again to FIG. 8C, at step 122, the chroma value of each interpolated patch is calculated based on the specified colorant values. Any conventional technique for determining chroma values based on colorant values may be used. For example, the well-known classical Neugebauer equations may be used to calculate the chroma values for each interpolated patch. At step 124, chroma values are calculated for the multi-colorant patches printed in step 110, and illustrated in FIGS. 10E-10G, using the same technique used in step 122. Thus, for example, if the classical Neugebauer equations were used at step 122, they are also used at step 124. Next, at step 126, any deviation is determined between the chroma values calculated at step 124, and the chroma values of the multi-colorant patches measured at step 112. At step 128, the chroma values of the interpolated patches calculated at step 122 are corrected using any deviation determined at step 126. For example, well-known vector-corrected Neugebauer equations may be used to correct the calculated chroma values. FIGS. 10H-101 illustrate vector-corrected chroma values for each of the interpolated patches. Finally, at step 130, the interpolated patch having the minimum vector-corrected chroma value is identified. In the example illustrated in FIGS. 10H-101, the minimum vector-corrected chroma value is 0.13, corresponding to an interpolated patch having CMY values (47, 56, 85). Thus, the optimal cyan, magenta and yellow colorant values that provide the best gray-balance are: C_(OPT)=47%, M_(OPT)=56% and Y_(OPT)=85%.

Referring again to FIG. 8C, at step 126, the optimal limit K_(OPT) is determined for the K colorant, which is equal to a weighted average of the minimum and maximum limits. An exemplary optimal limit for K may be expressed as:

K _(OPT)=(α×K _(MIN))+(β×K _(MAX))  (7)

where α and β are weighting factors. Exemplary values for the weighting factors are α=0.7 and β=0.3. Thus, from the exemplary minimum and maximum limits K_(MIN)=42% and K_(MAX)=69%, respectively, K_(OPT)=50%. Persons of ordinary skill in the art will understand that other values may be used for α and β.

Linearization Table Calculation

Referring now to FIGS. 3 and 11, exemplary methods and apparatus in accordance with this invention are described for calculating colorant linearization tables. In particular, FIG. 11 illustrates an exemplary method for calculating a colorant linearization table for a single colorant (e.g., Y). Persons of ordinary skill in the art will understand that the method may be repeated for each colorant used by proofing printer 46. Beginning at step 150, proofing printer 46 is used to print test pattern 82 including test patches P on output page 74. An exemplary output page 74 c including exemplary test pattern 82 c is illustrated in FIG. 12. Test pattern 82 c includes two strips of test patches, with each strip including fifteen test patches P. Persons of ordinary skill in the art will understand that test pattern 82 c may include more or less than two strips, and each strip may include more or less than fifteen test patches P. Each test patch P is comprised of a corresponding specified percentage of the colorant for which the linearization table is being created (Y in this example).

FIG. 13 illustrates exemplary colorant values (in percent) for yellow test patches P (patches are identified in each table by row (A-B) and column number (1-15)). The exemplary values provide test patches P having scales of single-colorant values. Persons of ordinary skill in the art will understand that other specific colorant values also may be used for test patches P. Referring again to FIG. 11, at step 152, calorimetric values are determined for each test patch P printed in step 150. For example, measurement device 76 may be used to determine XYZ data for each test patch P on output page 74 c. FIG. 13 illustrates exemplary measured XYZ data for test patches P.

Referring again to FIG. 11, at step 154, tonal response values are calculated from the measured calorimetric values from step 152. As previously mentioned, X-values are inversely proportional to amounts of cyan colorant, Y-values are inversely proportional to amounts of magenta colorant, and Z-values are inversely proportional to amounts of yellow colorant. In addition, Y-values are inversely proportional to amounts of black colorant. As a result, the tonal response of a cyan colorant patch P of percent i may be expressed as:

$\begin{matrix} {{{Tonal}\mspace{14mu} {Value}} = \frac{\left( {1 - \frac{X_{i}}{X_{W}}} \right)}{\left( {1 - \frac{X_{100}}{X_{W}}} \right)}} & \left( {8a} \right) \end{matrix}$

where X_(i) is the X-value for patch P, X_(W) is the X-value of paper white, and X₁₀₀ is the X-value of solid cyan (i.e., C=100%).

Similarly, the tonal response of a magenta colorant patch P of percent i may be expressed as:

$\begin{matrix} {{{Tonal}\mspace{14mu} {Value}} = \frac{\left( {1 - \frac{Y_{i}}{Y_{W}}} \right)}{\left( {1 - \frac{Y_{100}}{Y_{W}}} \right)}} & \left( {8b} \right) \end{matrix}$

where Y_(i) is the Y-value for patch P, Y_(W) is the Y-value of paper white, and Y₁₀₀ is the Y-value of solid magenta (i.e., M=100%). Likewise, the tonal response of a yellow colorant patch P of percent i may be expressed as:

$\begin{matrix} {{{Tonal}\mspace{14mu} {Value}} = \frac{\left( {1 - \frac{Z_{i}}{Z_{W}}} \right)}{\left( {1 - \frac{Z_{100}}{Z_{W}}} \right)}} & \left( {8c} \right) \end{matrix}$

where Z_(i) is the Z-value for patch P, Z_(W) is the Z-value of paper white, and Z₁₀₀ is the Z-value of solid yellow (i.e., Y=100%). Similarly, the tonal response of a black colorant patch P of percent i may be expressed as:

$\begin{matrix} {{{Tonal}\mspace{14mu} {Value}} = \frac{\left( {1 - \frac{Y_{i}}{Y_{W}}} \right)}{\left( {1 - \frac{Y_{100}}{Y_{W}}} \right)}} & \left( {8d} \right) \end{matrix}$

where Y_(i) is the Y-value for patch P, Y_(W) is the Y-value of paper white, and Y₁₀₀ is the Y-value of solid black (i.e., K=100%). FIG. 13 illustrates exemplary tonal values calculated based on measured XYZ data for yellow test patches P. FIG. 14 illustrates a graph of the calculated tonal values versus input values, with curve 170 representing the calculated tonal values of FIG. 13.

Referring again to FIGS. 3 and 11, at step 156, XYZ data are retrieved from backwards transform 32 b of press profile 32 for each colorant for 0%, 100%, and a midtone value (e.g., 40%, 50% or other similar midtone value). Table 2 illustrates exemplary XYZ values extracted from press profile 32:

TABLE 2 C M Y K X Y Z 100 0 0 0 15.22 22.09 55.16 0 100 0 0 32.36 16.37 15.04 0 0 100 0 68.41 73.38 5.98 0 0 0 100 2.22 2.36 1.94 40 0 0 0 52.71 58.79 70.22 0 40 0 0 60.90 53.95 51.39 0 0 40 0 76.98 81.94 43.15 0 0 0 40 39.55 41.03 37.38 0 0 0 0 85.67 88.32 78.70

Next, at step 158, the tonal response of press 34 is calculated using the XYZ data retrieved from press profile 32 in step 156, and equations 8(a)-8(d). Using the exemplary values in Table 2, the tonal response for yellow is:

TABLE 3 Input Value (%) Tonal Value (%) 0 0 40 49 100 100

At step 160, the tonal response of press 34 is calculated for the entire range of input values from 0-100% based on the tonal values calculated in step 158. For example, a spline function may be used to calculate tonal response values for the entire range of input values using the three data points in Table 3. Curve 172 in FIG. 14 represents the tonal response values calculated in step 160. Referring again to FIG. 11, at step 162, linearization tables are created by mapping press input values to equivalent printer input values using the tonal response data. For example, using tonal response curves 170 and 172, input values of press 34 are mapped to corresponding input values for proofing printer 46 that have equivalent tonal responses. The process of FIG. 11 may be repeated for each colorant used by proofing printer 46.

Distribution of Multi-Hue Colorants

If proofing printer 46 uses light and normal hues of a colorant (e.g., cyan), the printer may use only light cyan over a first range of input values (e.g., 0-100%), and may use a combination of light and normal cyan over a second range of input values (e.g., 40-100%). Referring now to FIGS. 3 and 15, exemplary methods and apparatus in accordance with this invention are described for determining distribution functions of multi-hue colorants for such systems. Beginning at step 180, processor 78 determines the lower limit L1 of the second range of input values (e.g., 40%), and the amount A1 of light colorant at 100% (e.g., 5% light colorant). For example, processor 78 may prompt a user to provide these two values.

Next, at step 182, proofing printer 46 is used to print test pattern 82 including test patches P having scales of light colorant only, scales of normal colorant only, and combinations of light and normal colorant. Each test patch P is comprised of a corresponding specified percentage of light and normal colorants. Next, at step 184, calorimetric values are determined for each test patch P printed in step 182. For example, measurement device 76 may be used to determine XYZ data for each test patch P. A step 186, tonal response values are calculated from the measured colorimetric values from step 184. As previously mentioned, for cyan, tonal values are calculated using equation 8(a), for magenta, tonal values are calculated using equation 8(b), for yellow, tonal values are calculated using equation 8(c), and for black, tonal values are calculated using equation 8(d).

Next, at step 188, the tonal response of press 34 is calculated as in steps 156-160 of FIG. 11, thereby providing target tone values for the multi-hue colorant from 0-100%. At step 190, the light distribution function is calculated from 0% to L1 as in step 162 of FIG. 11. Referring again to FIG. 15, at step 192, the light colorant distribution between L1 and 100% is determined using any suitable curve fitting techniques based on the slope of the light colorant distribution curve at L1, the slope of the light colorant distribution curve at 100% (0).

Next, at step 194, a two-dimensional table is created that provides the calculated tone value for each combination of light and normal colorant calculated at step 186. Finally, at step 196, the normal colorant distribution from L1 to 100% is calculated by determining from the two dimensional table the amount of normal that when added to the light value gives the target tone value. Persons of ordinary skill in the art will understand that the light and normal colorant distribution curves may be smoothed using any conventional smoothing algorithm.

The foregoing merely illustrates the principles of this invention, and various modifications can be made by persons of ordinary skill in the art without departing from the scope and spirit of this invention. 

1. A method for calibrating a digital color imaging device to a printing press by determining a plurality of distribution functions for a plurality of multi-hue colorants, the method comprising: printing a test pattern comprising a plurality of patches, each patch comprising a plurality of colorimetric values; measuring a calorimetric value for each patch; calculating a tonal response from said measured colorimetric values; calculating a tonal response of said printing press; calculating a first light distribution; calculating a second light distribution; creating a table of light-normal tonal responses; and calculating a normal distribution.
 2. The method of claim 1, wherein said tonal response of said printing press is from 0 to 100%.
 3. The method of claim 1, wherein said first light distribution is from 0% to a lower limit.
 4. The method of claim 1, wherein said second light distribution is from a lower limit to 100%.
 5. The method of claim 1, wherein said normal distribution is from a lower limit to 100%.
 6. The method of claim 1, wherein said step of measuring a colorimetric value is performed by a spectrophotometer.
 7. The method of claim 1, wherein said step of measuring a colorimetric value is performed by a spectrocolorimeter.
 8. The method of claim 1, wherein said test patches comprise scales of light colorant only, scales of normal colorant only, and combinations of light and normal colorant.
 9. The method of claim 1, wherein said colorant comprises cyan.
 10. The method of claim 9, wherein said tonal response value for cyan is calculated using the following equation: ${{Tonal}\mspace{14mu} {Value}} = \frac{\left( {1 - \frac{X_{i}}{X_{W}}} \right)}{\left( {1 - \frac{X_{100}}{X_{W}}} \right)}$ wherein X_(i) is an X-value for said patch, X_(W) is a X-value of a white paper, and X₁₀₀ is an X-value of solid cyan.
 11. The method of claim 1, wherein said colorant comprises magenta.
 12. The method of claim 12, wherein said tonal response for magenta is calculated using the following equation: ${{Tonal}\mspace{14mu} {Value}} = \frac{\left( {1 - \frac{Y_{i}}{Y_{W}}} \right)}{\left( {1 - \frac{Y_{100}}{Y_{W}}} \right)}$ wherein Y_(i) is a Y-value for said patch, Y_(W) is a Y-value of a white paper, and Y₁₀₀ is a Y-value of solid magenta.
 13. The method of claim 1, wherein said colorant comprises yellow.
 14. The method of claim 13, wherein said tonal response for yellow is calculated using the following equation: ${{Tonal}\mspace{14mu} {Value}} = \frac{\left( {1 - \frac{Z_{i}}{Z_{W}}} \right)}{\left( {1 - \frac{Z_{100}}{Z_{W}}} \right)}$ wherein Z_(i) is a Z-value for said patch, Z_(W) is a Z-value of a white paper, and Z₁₀₀ is a Z-value of solid yellow.
 15. The method of claim 1, wherein said colorant comprises black.
 16. The method of claim 15, wherein said tonal response for black is calculated using the following equation: ${{Tonal}\mspace{14mu} {Value}} = \frac{\left( {1 - \frac{Y_{i}}{Y_{W}}} \right)}{\left( {1 - \frac{Y_{100}}{Y_{W}}} \right)}$ wherein Y_(i) is a Y-value for said patch, Y_(W) is a Y-value of a white paper, and Y₁₀₀ is a Y-value of solid black.
 17. A system for calibrating a digital color imaging device to a printing press by determining a plurality of distribution functions for a plurality of multi-hue colorants comprising: a color output device for printing a test pattern comprising a plurality of patches, each patch comprising a plurality of colorimetric values; a measurement device for measuring a calorimetric value for each patch; a processor for calculating a tonal response from said measured colorimetric values; a processor for calculating a tonal response of said printing press from 0 to 100%; a processor for calculating a first light distribution from 0% to a lower limit; a processor for calculating a second light distribution from a lower limit to 100%; a processor for creating a table of light-normal tonal responses from a lower limit to 100%; and a processor for calculating a normal distribution.
 18. The system of claim 17, wherein said measurement device is a spectrophotometer.
 19. The system of claim 17, wherein said measurement device is a spectrocolorimeter.
 20. The system of claim 17, wherein said colorants comprise magenta, yellow, cyan, and black. 